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Abstract 

We show that R-parity violating decay of Wino dark matter of mass ~ 3 TeV 
can naturally account for the flux and spectral shape of the cosmic-ray electrons and 
positrons observed by the PAMELA and Fermi satellites. To provide a theoretical 
basis for the scenario, we also present a model that trilinear R-parity breaking ap- 
pears with a coefficient suppressed by powers of the gravitino mass, which naturally 
leads to the Wino lifetime of 0(10^^) sec. 



1 Introduction 



The lightest supersymmetry (SUSY) particle called as LSP in the SUSY standard model 
(SSM) is known as a good candidate of dark matter (DM) in the Universe if its mass 
is O(IOO) GeV — 0(1) TeV. The stability of the LSP can be guaranteed by assuming an 
exact R parity. However, the R parity may not necessarily be an exact symmetry. In 
fact, the LSP can be still DM as long as R-parity breaking terms are sufficiently small 
and the lifetime of the LSP is much longer than the present age of the Universe. If this 
is the case, the decay of DM may give rise to some excesses in cosmic rays. 

Much attention was recently attracted to the anomalies in cosmic-ray electron/positron 
fluxes observed by PAMELA [1] and ATIC [2], and decaying LSP scenarios were exten- 
sively discussed in this context PiUgE Most of the proposals assume bilinear R-parity 
breaking terms such as LHu-, since they are the lowest dimensional R-parity breaking 
operators which most likely dominate R-parity breaking effects at low energies. However, 
the bilinear R-parity breaking terms induce DM decays into quarks and hence produce 
too many antiprotons in cosmic rays. Therefore, it is much safer to consider the next 
lowest dimensional operators such as eLL [5l [TJ El [9] . However, it seems very unnatural 
to consider the second lowest dimensional operators, suppressing the lowest dimensional 
bilinear operators. In addition, even if the trilinear term dominates over the bilinear one, 
the magnitude of the R-parity breaking should be extremely small, especially in the case 
of neutralino LSP. No rigorous explanation for the smallness was known, and the size of 
the R-parity violation was treated as a free parameter. 

Very recently, the Fermi collaboration has released data on the electron/positron fluxes 
from 20 GeV up to ITeV [10]; the spectrum falls as E^^'^ without prominent spectral 
features, and it is in ag reement with the H.E.S.S. data at E ~ 1 TeV [II1II2]@ The index of 
the observed electron/positron spectrum is close to the high end of theoretically expected 
value. Moreover, if we combine the Fermi data for E'^lTeV and the H.E.S.S. data 
for E^l TeV, it looks that the spectrum becomes softer at energies above 1 TeV. From 
this observation, we regard the relatively hard (and almost featureless) electron/positron 

^ See also Ref. [6] for other decaying DM models. 

^ The Fermi result is not consistent with the excess around E ~ 600 GeV reported by ATIC. In this 
letter we adopt the Fermi result and do not consider the ATIC data. 
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spectrum below ITeV reported by Fermi as an excess with respect to the background. 
If the electron/positron spectrum observed by Fermi (as well as the positron fraction 
observed by PAMELA) is to be explained by the DM decay , the mass of DM must be 
about a few TeV. The suggested mass scale is intriguingly close to the mass of a Wino 
LSP required for the thermal relic to explain the observed DM abundance. 

In this letter we present a model that trilinear R-parity breaking appears with a 
coefficient suppressed by powers of the gravitino mass. Interestingly, the size of the R- 
parity breaking is naturally small, and moreover, it leads to the required size of the 
R-parity breaking suggested by observation, in the case of the Wino LSP of a mass 
~ 3 TeV. Based on this model, we show that the Fermi as well as PAMELA data can be 
simultaneously explained by the Wino LSP decaying through the trilinear term eLL. 

2 A model of R-parity breaking 

It is known that the superpotential possesses a constant term Co to cancel the positive 
energy density induced by SUSY breaking. The constant term is equal to the gravitino 
mass, Co = "023/2, in the Planck unit in which the Planck mass Mp ~ 2.4 x 10^^ GcV is set 
to be unity. The constant term breaks the continuous U{1)r symmetry down to a discrete 
Z2 symmetry, which is nothing but the R parity. A crucial observation is that, if the R 
symmetry at high-energies is not the continuous U{1)r but a discrete Z2k+i {k = integer), 
the constant term Co results in the R-parity breaking. In this letter we take k — 2, since, 
as we will see later, it leads to the required magnitude of the R-parity breaking operators 
to account for the anomalous excess in cosmic-ray electrons and positrons. 

Let us take the R charges for all quark and lepton chiral multiplets to be 1 and those 
for the Higgs chiral multiplets and to be 0. (See Table 1.) Then, the following 
lepton-number violating trilinear term is allowed by the Z^r symmetry: 

W = Hijk{CQ)%LjLk, (1) 

where Kijk is a numerical coefficient, i, j, /c = 1, 2, 3 denote the generation, and summation 
over the SU(2) gauge indices is understood. The coupling Kijk must be antisymmetric 
under the exchange of the last two indices {j <-> k) due to the SU(2) gauge invariance. 
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Figure 1: One of the diagrams representing the Wino decay through the trihnear R-parity 
breaking ([1]). 

Note that the R charge of Co(= ^3/2) is 2. We may presume that the coefficient Hijk 
takes a larger value for the third (and second) generation. Therefore, in the following we 
focus on Kj23 with i = 1,2, 3, which is assumed to be unsuppressed with the other terms 
of different combination of flavors being suppressed. 

We assume that a Wino LSP accounts for the DM in the Universe. In order to explain 
the observed DM abundance, the Wino mass must be in the range between 2.7—3 TeV |13j . 
The Wino LSP is naturally realized in the anomaly-mediated SUSY breaking [14] . Then, 
using the relation between the Wino mass and the gravitino mass in the anomaly medi- 
ation, we find that the gravitino mass should be about 10^ TeV. In the presence of the 
R-parity breaking ([1]), the Wino LSP is no longer stable, and decays into neutrinos and 
charged leptons through the exchange of a virtual slepton. One of the decay diagrams 
is shown in Fig. [TJ The decay rate of the Wino LSP through the interaction ([T]) with 
(j, k) = (2, 3) is given by [I5] 

V ^^ ^ , tH' ^ J \ J \ ^2s\ ^^lo^TeV/ V3 TeV/ V5TeV/ ' 

where denotes the Wino mass, and we have assumed the common slepton mass, m^, 
for simplicity. As we will see in the next section, the lifetime ([2]) is close to what is needed 
to explain the cosmic-ray observation. 

Several comments are in order. First, there could be other trilinear terms including 
quark multiplets, but they may have only negligible effects on the decay processes if 
squarks are substantially heavier than sleptons. Second, the bilinear term is also allowed 
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Table 1: The assignment of R- charge and Z^. 



by the Z^r symmetry and takes the following form: 

W = [C^fLH^. (3) 

Since the bilinear term is suppressed by an additional factor Cq compared to the trilinear 
term, the latter becomes much more effective than what is naively expected based on the 
dimensional grounds. It is not trivial, though, if the trilinear term dominates over the 
bilinear term as the decay processes of the LSP. We will come back to this issue in Sec. HJ 
and for the moment, we simply assume that the bilinear term does not have sizable effects 
on the decay. 

3 Cosmic-ray signal from LSP decay 

Let us discuss the cosmic ray signals from the LSP decay. Its decay pattern depends 
on R-breaking structure and the SSM mass spectrum. For a demonstration, we consider 
the case that the eiL2L^ {i = 1,2,3) term dominates the R-breaking and assume that 
BF(DM — > T^u^ef) = BF(DM — » u^fi'^ef) = 0.5. We have used the constant matrix 
element in the three-body phase space, for simplicity. The electron and positron energy 
spectrum is estimated with the program PYTHIA |16] . For the propagation of the cosmic 
ray in the Galaxy, we adopt the same set-up in Ref. [IT], based on Refs. fiSl fT9] . As for 
the electron and positron background, we have used the estimation given in Refs. [20l[2T] . 
with a normalization factor fcbg = 0.68. In Fig. [2l we show the positron fraction and the 
electron and positron total flux. We set that ttidm = 3 TeV and the lifetime is 9 x 10^^ 
sec. As can be seen from Fig. |2l the cosmic-ray signal in the present model can nicely 
fit the PAMELA data for i = 1, 2, 3. On the other hand, the prediction in case of i = 1 
fails to explain the Fermi data due to the presence of hard electrons produced by the LSP 
decay, leading to the bump at 1 TeV. The prediction in case of i = 2 gives a very good fit 
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Figure 2: Cosmic ray signals in the present model, (a): positron fraction with ex- 
perimental data [U l22l [23]. (b): positron and electron fluxes with experimental 
data [ini [H], [121 [21 [21] . The yellow zone shows a systematic error and the dashed line 
shows the background flux. I, II and III represent the cases that eiL2L^, 62-^2-^3 and 
631/2-^3 dominate the R-breaking, respectively. 



to the Fermi data. Considering uncertainties in the background estimation as well as the 
SSM mass spectrum, however, the case of i = 3 may be able to give an equally good flt. 

4 Discussion and conclusions 

Let us discuss several issues in the model presented in Sec. [21 and possible solutions. First 
of all, the bilinear term like ([3|) could induce additional decay processes into lepton and 
Higgs, which may result in too many antiprotons. Note that the effective mixing angle 
induced by the bilinear term is given by the ratio of Cq to the Higgs mass or slepton mass. 
Therefore, if m^i2 is the same orders of magnitude of the Higgs mass or the slepton mass, 
the bilinear term will be as important as the trilinear term. 

To make the matter worse, the bilinear term will be enhanced in the presence of right- 
handed neutrinos A^. This is because, if we assign the R-charge 1 to as in Table 1, the 
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following interactions are allowed by the Z^r symmetry. 

W = C^N + ^MNN + y.NLHu, (4) 

where M denotes the Majorana right-handed neutrino mass, y^, is the neutrino Yukawa 
coupling, and we have suppressed the flavor indices. Since the first term is generically 
present, the right-handed neutrino will develop a non- vanishing expectation value of 
0{Cq/M). Then, the neutrino Yukawa coupling induces the bilinear term y (N) LHu ~ 
{C^/M)LHu, which is enhanced by 1/M compared to ([3]). Thus, the presence of the 
right-handed neutrino of mass lighter than the Planck mass leads to the enhancement of 
the bilinear term. 

Those problems can be easily solved by introducing a symmetry, under which the 
SSM particles and M are charged by an unit charge. (See Table 1.) Namely, the Z3 
symmetry is broken by M <^ 1. Then the bilinear term is suppressed by M, while the 
first term in (j4]) by M^. Therefore, the trilinear term ([1]) will be the dominant source of 
the LSP decay. 

In this letter we have presented the R-parity breaking model in which the trilinear 
term becomes important and appears with a coefficient proportional to powers of the 
gravitino mass. We have also shown that the Wino LSP decaying through the lepton- 
number violating trilinear term ([T]) can account for the PAMELA and Fermi data. 

There are several non-trivial coincidence in our scenario. First of all, the change 
in the power index of the electron/positron spectrum suggests the DM of mass a few 
TeV. This is intriguingly close to the mass of the thermal relic Wino DM. Assuming 
the anomaly mediation, the gravitino mass is determined to be about 10^ TeV. In our 
model on the R-parity violation, the lifetime of the Wino LSP is determined by some 
combination of the Wino mass and the gravitino mass. Substituting the above values 
for the Wino and gravitino masses, we have obtained the lifetime which is surprisingly 
close to what is needed to account for the cosmic-ray anomalies. Further study of the 
prediction of other cosmic rays such as gamma-rays, antiprotons and neutrinojfl, as well 
as future observational data, will enable us to tell whether those are just coincidence or 

Our model is consistent with the current observational bound on the neutrino production from 
decaying DM [25] . 
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Figure 3: Predicted signals of diffuse gamma-ray flux shown together with the EGRET 
data [271 EH]. 

may reflect from the characteristics of DM and the underlying physics beyond the SM. 

We can see from Fig. [2] that the fit to the Fermi data becomes better if the decay 
products include the charged leptons in the second and third generations, namely, muons 
and taus. In particular, if the decay product is dominated by the muon (as in the case II 
of Fig. [2]), the fit looks pretty goodjj It is actually possible to give a equally nice fit to the 
Fermi data, if we properly combine the contributions from the first and third generations. 
This is indeed the case if we consider unsuppressed ^313. Such flavor dependence may be 
probed by studying the diffuse gamma-ray in detail. For reference we show in Fig. [3] the 
predicted diffuse gamma-ray signals for the cases I, II and III. Including more taus in the 
final states generically lead to larger signal in the diffuse gamma-ray. We may be able to 
untangle the flavor dependence by making use of the different predictions. 
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